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Abstract 

We diagonalize the double-row transfer matrix of the SU{N) vertex model for certain 
classes of non-diagonal boundary conditions. We derive explicit expressions for the corre- 
sponding eigenvectors and eigenvalues by means of the algebraic Bethe ansatz approach. 
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The study of integrable models of statistical mechanics with arbitrary boundary conditions 
has gained a tremendous impulse after the work by Sklyanin jl]. This author has been able 
to generalize the quantum inverse scattering method j2] to include the important case of open 
boundaries. It turns out that an exactly solved spin chain with open boundary condition can 
be obtained through the following double-row transfer matrix 

T(A) = Tr^ [k^2\X)TAX)k':^\X)Tj\-X)] (1) 

where Ta{X) = Cj[l{X)C_al^i{X) . . .£_4i(A) is the standard monodromy matrix that generates 
the corresponding closed spin chain with L sites |2^. 

We recall that the symbol A denotes a A^- dimensional auxiliary space and A parameterizes 
the integrable manifold. The operator C_Aj{X) represents the bulk weights of the corresponding 
vertex model whose transfer matrix commutes with the spin chain Hamiltonian. The N x N 
matrices K^^\X) describe the interactions at the right and left ends of the open chain. One of 
the simplest integrable system is the fundamental SU {N) vertex model Pj whose Boltzmann 
weights are given by 

N N N 

£^,(A) = a(A) E e ® ^1 + KA) E eli^ ^ef,+ j: e^^J e^] (2) 

a = l a,,/3 = l a, 13 = 1 

where a(A) = A + 1, 6(A) = A and e\J^ are the usual Weyl matrices acting on the space V. 
A quite general class of open boundary conditions for this vertex model is represented by the 
following K-matrices jU IH] 



)(A) = Mif )dJ)(A) M^f> ^ D]^\X) = E e^^\X)e^^^ (3) 



r(±) 



N 



The elements of the Nx N matrices are arbitrary c-numbers and the functions ^^''^^(A) 
are given by 



U--A a = p+l,...,N l^+ + f + A a = p+l,...,N 



1 



where p is an integer with values on the interval 1 < p < N and ^± are free-parameters. Here we 
emphasize that each X-matrix ()3I4|) has altogether 2N — 1 arbitrary parameters characterizing 
the interactions at the appropriate boundary. 

The diagonalization of the transfer matrix ()1I2|) for general non-diagonal X-matrices fl3l4|) 
is a tantalizing problem due to the difficulty of finding a suitable reference state to perform 
a Bethe ansatz analysis. However, progress on this matter has recently been done in the 
literature, most of it concentrated on the eight [Zj and six jSl El d]] vertex models. The 
Belavin model is to our knowledge the only multistate vertex system investigated so far 
with non-diagonal open boundaries jTT]. Though its bulk weights are known to reduce in 
the isotropic limit to those of the SU{N) vertex model the same does not occur for the 
boundary K-matrices. In fact, the elliptic i^'-matrices associated to the Zjy Belavin model J2] 
have fewer free-parameters, which totals + 1, as compared to that contained in the isotropic 
X-matrices ()3f4|l . Therefore, for > 3 the SU{N) vertex model with open boundaries is 
indeed a genuine integrable system that deserves to be studied independently. We suspect 
that this situation extends to many isotropic integrable vertex models based on higher rank 
symmetries. 

The purpose of this work is to show that the diagonalization of the transfer matrix (|H2j) 
of the SU{N) vertex model in the case can be mapped on a similar eigen- 

value problem with the diagonal boundaries D^\X). This constraint does not imply that the 
right and left fT-matrices are the same because the parameters ^± are still unrelated. This 
observation not only allows us to solve the eigenvalue problem for 2N independent boundary 
parameters but also makes it possible the relation between eigenvectors of seemly different 



open boundaries. In order to see that we insert the terms ■* 
the double- row transfer matrix (Q), permitting us to rewrite it as 



all over the trace of 



T(A) = Tr^ [D^J:\X)fAX)D[,\X)f^\-X)] (5) 
where the new monodromy 74(A) = >C^l(A)£^l_i(A) . . .£_4i(A) whose gauge transformed £- 
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operators are given by 



M 



(-) 

A 



(-) 



(6) 



Further progress is made by reversing the gauge transformation (jH)) with the help of the 
following quantum space transformation 

\Jt^Ia,[\)\J, = CA,i\) (7) 



where Ui = Id( 



IdOM 



(-) 



)Id, 



)Id and Id is the N x N identity matrix. This remarkable 



-th 



property [12] can now be used to define a new double-row operator T{X) 

TW = I[Ui'nX)f[U, (8) 

By using the canonical transformation (jH)) together with the assumed constraint = 
between the right and the left matrices we find that T(A) becomes 



T(A) = Tr^ DW(A)r^(A)D^-^(A)r^-i(-A) 



(-), 



(9) 



which is precisely the double-row transfer matrix of the SU{N) vertex model with diagonal 
X-matrices D^^\X). 

As a consequence of that the operators T{X) and T{X) share the same eigenvalues and 
furthermore if ■0^ is an eigenstate of T(A) then the corresponding eigenvector 1-0) of T{X) 

L 

is Y\_ Uj if^. From now on our main task consists therefore in diagonalizing the double-row 

transfer matrix T'(A). In this case the associated open boundaries D^'*(A) are diagonal and 
such eigenvalue problem can be tackled by standard nested Bethe ansatz approach. By now 
this procedure has been well explained in the literature, see for instance refs. ^JE]; and here 
we shall restrict ourselves in presenting only the essential steps of the solution. We first note 
that diagonal boundaries permit us to use as pseudovacuum the usual ferromagnetic state 



n®io),' io), 







(10) 



N 



The next step is to write convenient commutation rules for the elements of the double 
transition operator T^(A) = Tj[{X) Dj^\X)T_^^ {— X) which satisfies the following quadratic 
relation 



'11] 



We proceed by looking for a representation of the operator T^(A) that is capable to dis- 
tinguish potential creation and annihilation fields over the state |0). Previous experience with 
nested Bethe ansatz diagonalization of the SU (N) vertex model suggests us the form 



/ 



^(A) Si(A) 
Ci{X) Dn(A) 



Bn-i{X) 
DiN-iiX) 



\ 



(12) 



\ C;v-i(A) D^-ii(A) ■■■ Dr,-m-iiX) J 

From Eqs. ()lipi2j) it follows that out of A^^ possible commutation rules there exists three 
families that are of great use, namely 
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A{X)B,{^i) 



^(^-^)^^^ + ^)5,(^)A(A) ^^^^^ 



— A) a{fi + A) 

\ B,{X)D,,{^) 
a{X + fi) 



a(2/i) b{ii — A) 



i?,(A)A(^) 



(13) 



f«(A + /x + l) f«(A-/i) 



id 



b{X + 12 + l)b{X - fl) 



'^Bd{fi)D,g{X) 



f«(2A + l) 



id 



a{2X)b{X - fj,) 



+ 



6(2/i) 



fW(2A+ 1) 



id 



kj 



B,{X)B,{^) 



a{2fi) a{2X)a{X + fi) 
f«(A-/i) 



BdiX)Aifi) 



Bd{X)Dek{fi) 



(14) 



a(A — /i) 



Bc{fi)B4X) 



(15) 



where 



r^^\X) '' are the matrix elements of the SU{N — 1) /^-operator ^ (j2)) and the new Dij{X) 



kl 



operators are conveniently defined by the relation 

A,(A) = A,(A)- 



a(2A) 



A{X) 



(16) 



^Here we recall that we have used the convention (A) — |^f''-'(A) 
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abed 



Yet another important element is the action of the fields A{X), Dij{X) and Cj(A) on the 
reference state |0). This can be computed by using the triangularity property of Cj,j{X) upon 
\0)j and they are given by 



A{X) |0) 
A,(A) |0) 
Q{X) |0) 







2L 



4;! (A) - 



a(2A) 



% |0) 



;i7) 



We have now gathered the basic ingredients to perform an algebraic Bethe ansatz analysis. 
We suppose that further eigenstates 'i/'m) of T{X) can be put into the following structure 

= B,,iX^^^) . . .5a„,(AS)^--- |0) (18) 

where the indices aj run over A^ — 1 possible values and the rapidities {A^-^^} will be determined 
by solving the eigenvalue equation 



N-l 



a(2A) 



a=l 



N-l 



MX) V'ni) + E e^a+lWDaa{\) ^n, ) = A(A) Vni ) (19) 



By carrying on the fields ^(A) and Dii{X) over the multiparticle state psp we generate 
terms that are proportional to ipm'^ and those that are not, denominated unwanted terms. 
The first ones contribute to the eigenvalue A(A) and are obtained by keeping only the first 
terms of the commutation rules ()13II15|) and by requiring that the coefficients jF""i --'^i are 



eigenstates of an auxiliary double-row operator T*^^)(A, {A|-^^}) given by 



f«(A,{A«}) = Tr^,.) 



such that 

^i;i(A,{A«}) = f%,^iX + X? + l)...f%,^JX + Xil+l) 
fjilMAXf}) = r«,,„/A-A«J...f«,,(A-A«) 

where A^^^ G C^^^ and the associated iT-matrices are 

e[-\X) 



(20) 



Af-l 



Af-1 



D'lltdW = E ei(A) elf' D%\X) = E 



a=l 



a=l 



a(2A) 



(AW) 



(21) 



(22) 



The diagonalization of the inhomogeneous operator T^^\X, {X^p}) is implemented by a 
second Bethe ansatz which is once again parameterized by a new set of rapidities Ap"* . . . A^^^'*- 
By keeping on going this procedure we are able to relate the eigenvalues A*^'^(A, {A^'"*}) of the 
transfer matrix r('^(A, {A^''}) at the nearest neighbor steps / and / + 1. Since the commutation 
relations for the elements of the corresponding double transition operator is similar to that 
exhibited in Eqs. p3lll5p it is not difficult to derive the following recursive relation 

A'«(A, {Af}) - ft «(A . AS" . 0«(A - A?') H Z ^ g t g! I g 

1 



i=l 

ni "i+i 



+ UKX + xf' + mX-Xf^) n , , , , ^^ A^^+^)(A,{Ar)}) 



(23) 



where the functions Q^''\X) are given by 



I (A+'haJ'+M + A) (f - P + e+ - A) / = 0,...,p-l 



ry(e--A-p)(f + e+ + A) /=p,...,iV-l 



(A+i)(A+i±i) 

For sake of consistency with our original eigenvalue problem we set A^"'* = for j = 
1, . . . ,no = L. By the same token the unwanted terms generated in the eigenvalue problem of 
the double-row operator T*^'''(A, {A^^}) are cancelled out provided that the variables {Xf^^^} 
satisfy the nested Bethe ansatz equations 

"f^^ a{4^ + Af-^^ + / - 1) a(Ai'^ - xf''^) Q^^-'Kx'P) &(2A? + / - 1) ^ 
ii KXi^ + Xf-'^ + 1-1) biXl^ - Xf^) gW(Ai')) a(2A? + /) 

L{ KAi'^ + Af + / - 1) a(Af - A?) M KAf^) - A?) a{x!^-'^ + A« + /) 

(25) 

Explicit results are now obtained by iterating Eqs. p3p25j) beginning at / = until we reach 
the step I = N — 2. At such final step one has to diagonalize an inhomogeneous six vertex 
model with open boundaries by adapting previous results obtained by Sklyanin p. Putting 
together all that and by making the convenient displacements A^''' Xf'' — | we find that the 
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final result for the eigenvalues A (A), up to a normafization factor of value (1 — A^)^, are given 
by the expression 

A(A) = Q(°)(A) [a(A)]^^ fl ^ " D + ^i"^ " i) 



.=1 (A - Af) + i) (A + + I 



' ^ " h ^ ^ ' l\ (A - A»> + i) (A + Af + i) l\ (A - Ar " + 1±1) (A + Ar" + •-^) 

while the Bethe ansatz roots {A-^-*, . . . , A-^^"^'} satisfy the following system of non-hnear 
equations 



(Afc ^ - |) 



2L 

(1) 



e«(Ai^^) = n 



X 





-A« 


+ 1) (Ai" 


+ A« 


+ 1) 




"t 


- 1) (At" 


+ A« 


-1) 


(Af 


~ Afc 


+ i) (Af ' 


+ At" 


-^) 


(Af 


-Al" 


- ^) (Af ' 


+ Af' 





112 



(Af-Af' + ^XAf + Ar' + i) ^,,,,.,,!, ^ ^ (AW-Af + lXAf+Af + 1) 
(Ai« - a!'-" - i) (A« + Al'-'> - i) ' ,U (Ai« - Af - 1) (Ai" + Af - 1) 



X 



y (Af^^ - Af) + \) (Af^) + Aj!) - i) 

j=l (^Aj - A^, - 2 j l,Aj- + A^ + 2J 

1 = 2,. ..,N -2 (28) 



n l^fc ~ +2^ l^A: ^ \ + 2^ Q(iV-l) / A (A^-l) A 

..(7V-1) ,(iV-2) 1. , (^-1) (JV-2) iv'-' \^-k ) 

i=l I'^A; ~ ~ 2) l^fc + ^« ~ 2'' 

^ (Af - Af + 1) (Af + Af + 1) 
jy^fe (Afe - Aj- - Ij (Afe + Aj- - Ij 



where the function ©('^(A) is given by 

e(')(A) = <^ (A+|-«-) * - ^ (30) 

I 1 otherwise 

We would like to close this letter with the following remarks. The same strategy described 
above works when one of the boundaries is purely free and the other stays arbitrary with 
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2N — 1 independent parameters. This is for instance the case of i^^TlA) = Id and K'x^X) 
an arbitrary /C-matrix ()H|4j) . It turns out that the corresponding Bethe ansatz results for this 
choice of boundaries are obtained from Eqs. 1)2611^0)) by taking their ^ oo hmit. We note 
that similar result has recently been reported in ref. however on the basis of the analytical 
Bethe ansatz method and solely for case of diagonal boundaries. 

Although we have concentrated our attention on > 3 SU (N) models, similar idea is also 
applicable with success to the SU (2) XXX spin chain. In this special case it is more convenient 

all over the double-row transfer matrix (|ll2j) and after 
reversing the transformed /^-operators we impose that 
matrix ^. This allows us to carry out the algebraic Bethe ansatz for a single constraint between 
the six possible boundary parameters, reproducing what has been found earlier for the XXZ 
chain [0]. However, this method has the clear advantage of relating the eigenfunctions of the 
XXX chain having all possible boundary terms with five free parameters to that with only one 
off-diagonal and suitable diagonal boundary terms. We hope that this relationship could be 
of utility in physical applications such as in the study of the scaling behaviour of symmetric 
exclusion processes [T7] . 

The possibility of undoing gauge transformed ^-operators appears to be a general property 
of isotropic vertex models [13j. Therefore, we expect that the method devised in this work will 
be useful in the solution of the eigenspectrum of a variety of isotropic systems with certain non- 
diagonal open boundaries. Interesting examples would be the case of soliton non-preserving 
boundaries [inillH] for the conjugated representation of the SU (N) as well as general boundaries 
for vertex models invariant by the 0{N) and Sp{2N) Lie algebras. 

^This construction clearly leads us to fewer constrained boundary parameters than the condition = 
M_^~^ but when > 3 a Bethe ansatz analysis has eluded us so far. 



to start by inserting M^/' M^^^ 



W^A is a triangular 
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